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I. INTRODUCTION 

An adequate description for low-energy phenomena of fermions is provided by the mean- 
field approximation in which particles move under the mean-field Hamiltonian; relaxation 
effects are not taken into account in this approximation [1]. The relaxation effects due to the 
reservoir, which are important to the irreversibility [2], have been discussed in the literature 
[1-7]. It was shown that loss and gain should be considered for such effects and that the loss 
comes from the imaginary part of the self energy [4,8] in the quantum master equation[6]. 

To include the self energy, it is reasonable to extend the mean-field Hamiltonian to 
include nonhermitian operators. It will be shown in section II that the master equation for 
fermions can be derived from such an extension, if we assume that the equation is symmetric 
with respect to particles and holes. The corresponding equation for bosons is obtained in 
section III. The nonlinear master equation is obtained in section IV. To compare with other 
stochastic equations, it is shown in section V that the derived equations reduce to the 
Markoff master equation [3-5,7] in the low-density limit, and to the quasiclassical master 
equation [9- 11] for homogeneous systems. Conclusions are made in section VI. 

II. NONHERMITIAN OPERATORS FOR FERMIONS 

In the mean-field approximation, the system of identical particles is described by a one- 
particle density matrix pp{t) whose trace equals the number of particles. [1] To extend such 
an approximation, we consider the quantum relaxation effects. For fermions, it will be shown 
in this section that the equation, including the relaxation term in addition to the Liouville 
flow, should be of the general form, 

^^ppit) = -^[Hit),ppit)] + {ppit)^Mt)}-{i-Ppii)^Mt)}, (1) 

where Ap{t) and Ap{t) are hermitian operators arising from the relaxation effect, H{t) is the 
hermitian operator for the Liouville flow, and / is the identity operator. Here we take h = 1, 
and notations [A, B] and {A, B} denote, respectively, the commutator and anticommutator 



of operators A and B. As shown in section IV, from Eq. (1) we can intuitively obtain the 
relaxation term, which is formally derived in Appendix A from the second quantization. 
In the mean-field approximation, pp{t) is governed by 

Ppit) = U{t,t')p,{t')U\t,t'), (2) 

with the time-evolution operator U{t,t') generated by the mean-field Hamiltonian H{t), 
which is a hermitian operator. To include the self energy with imaginary part, it seems 
reasonable to extend H{t) to a Hamiltonian 7i(t) with a nonhermitian part, 

n{t) = Hit) + iA{t) (3) 

where A[t) is a hermitian operator. Then Eq. (2) is equivalent to 

^p,(t) = \[H{t),p,{t)] + {pp(t), A(t)}. (4) 

The first two terms on the right-hand side of Eq. (1) have been included in Eq. (4) when 
we set A{t) = Ap{t). However, if we assume the equation for fermions to be symmetric with 
respect to particles and holes, Eq. (4) must be modified to become Eq. (1), as shown later. 
The last two terms in Eq. (1) are for the loss and gain factors of particles. 

Instead of Pp{t), fermions can be described as well by holes, i.e., the vacencies of particle 
orbitals [10]. First, we diagonalize Pp{t) in an orthonormal basis {IV'a(^))}: 

Ppit) = J2c,it)\Mt)){Mt)l (5) 

A 

where c\{t) denotes the number of particles in \iljx{t)). Define 

Pp{t) = / - pp{t) = ^(1 - c,mMt)){Mt)\- (6) 

A 

For any normalized state \(p), 

(0Mt)|0) = l-(0|p,(t)|0), (7) 

which is just the number of holes in that state, and therefore it is reasonable to use Pp{t) to 
represent holes. In fact, we can obtain the density matrix of holes from Pp{t), as shown in 
Appendix B. In terms of Pp{t), Eq. (4) becomes 



^pp(t) = ^-[H{t),p^{t)] - {I-p^{t),A{t)}. (8) 

However, the form of the equation for Pp{t) is different from that for Pp{t) if A{t) ^ 0. 

It is known that in Eqs. (4) and (8), H{t) describes the Liouville flow. To see the physical 
meaning of A{t) in these two equations, let us consider the special case that H{t) = and 
A(t) = — -2 l^)(^l with 10) being a normalized state. Here 7(t) is real because A{t) is 
hermitian. From Eqs. (4) and (8) we obtain for particles and holes, respectively, 

^((/)|pp(t)|0) = -7W(0|PpWl0), (9) 

^(0|/>p(t)|0)=7(t)(l-(0lMt)|0)), (10) 

When j{t) > 0, Eq. (9) describes the loss of particles in |0). The loss of particles is just the 
gain of holes, and therefore Eq. (10) should describe the gain of holes in |0) if 7(t) > 0. 

It seems natural to set 7(t) < in Eq. (9) to describe the gain of particles in \(p). The 
resulting equation, however, is not of the same form as Eq. (10) with 7(t) > for the gain 
of holes. In fact, it seems natural to use 

|(0|PpWI0) = 7'(t)(i-(0Mt)|0)) (11) 

to describe the gain of particles, just as we describe the gain of holes using Eq. (10), where 
7'(t) is a positive real number. Similarly, if the loss of holes is described by Eq. (10) with 
7(t) < 0, the loss of particles should also be described by Eq. (11) with 7'(t) < 0. 

To correctly describe loss and gain of particles, we note that the loss rate of particles 
should be zero when |0) has no particle, while from Pauli effect, in fact, the gain rate of 
particles should be zero when the orbital is filled with particles. Therefore, Eq. (9) is only 
for the loss of particles with 7(t) > while the gain of particles is described by Eq. (11) with 
7'(t) > 0. Now the loss and gain rates are proportional to (0|pp(t)|0) and 1 — (0|pp(t)|0), 
respectively. In Appendix C it is shown that Eq. (4) cannot describe the gain of particles 
in 10) no matter what H{t) and A{t) are. 



Because Eqs. (9) and (10) are obtained from the second terms on the right-hand side 
of Eqs. (4) and (8), respectively, we can see why Pp{t) should be governed by Eq. (1), in 
which Ap{t) and Ap{t) are for the loss and gain of particles. Eq. (11) may be obtained from 
Eq. (1) with H{t) = Ap{t) = and Ap{t) = -2^|0)(0|, while Eq. (9) is obtained with 
H{t) = Ap{t) = and Ap{t) = -2^|0)(0|. We can see that Ap{t) is for the loss of holes 
just as Ap{t) for the loss of particles, by rewriting Eq. (11) with corresponding Pp{t), 

^(0|/^(t)|0) = -7'(t)(0|/^(t)|0). (12) 

Rewriting Eq. (1) with corresponding Pp{t), we arrive at 

-Ppit) = -^[Hit),pp{t)] + {pp(t),Ap(t)} - {I-p^(t),Apit)}, (13) 

which is symmetric to Eq. (1) for Pp{t). To include an antihertmitian part in addition to 
the mean-field Hamiltonian H{t), therefore, we need to include Ap{t) in addition to Ap(t). 
In section IV, the quantum master equation for fermions will be derived from Eq. (1). 
Although all Cx{t) in Eq. (5) are or 1 in the mean- field approximation, they could be 
fractions in the master equation [1,2]. 

III. THE RELAXATION TERM FOR BOSONS 

For a boson system described by the one-particle density matrix Pp{t), it is reasonable 
that the decay rate of particles is proportional to (0|pp(t)|0) for any normalized state |0). 
Unlike fermions, however, the gain rate for bosons may not vanish when (0|pp(t)|0) = 1, 
and hence we need to modify the equation for the gain of particles. By using the equation 

-p^(t) = -{I + ppit),Apit)}, (14) 

we have for the gain of particles in any arbitrary normalized ket \<p) 

^(0|p,(t)|0) = 7(t)(l + (0Mt)|0)), (15) 



where we set Ap{t) = — -2 l^)(^l with a real number 7(t) > 0. Now the rate for bosons to 
enter |0) is proportional to l + {(f)\pp{t)\(p). The enhancement of the gain rate with increasing 
{<f)\pp{t)\(f)) is consistent with the fact that bosons prefer states filled with many particles. 
Therefore, for bosons we shall consider the following equation 

-Ppit) = l[J/(t),p,(t)] + {ppit),Apit)} - {I + p^(t),Apit)} (16) 

to include the effects due to relaxation if Eq. (14) is suitable for the gain of bosons. As 
shown in the next section, from Eq. (16) we can intuitively obtain the relaxation term which 
is formally derived in Appendix A from the second quantization. Hence the modification for 
the gain of bosons is reasonable. 

IV. NONLINEAR QUANTUM MASTER EQUATION 

In the last two sections, we obtain Eqs. (1) and (16) so that the loss and gain factors 
are included in addition to the Liouville flow. To obtain the master equation, first let us 
consider a two-state system in contact with the reservoir and assume a time-independent 
Hamiltonian H{t) = Hq. Let \i) and |/) be normalized eigenkets of Hq. If the reservoir 
induces a transition so that particles in \i) tend to jump to |/), there should be a loss 
factor for \i) and a gain factor for |/). Therefore for Ap[t) and Ap{t) in Eqs. (1) and (16), 
respectively, we shall set 

A^Ht) = -^\f){fl (17) 



AHHt) = -^-^m^\, (18) 
where 71 (t) and 72 (t) are positive real numbers. So we have 

|(/|PpWI/)=7iW(l±(/|Pp(t)|/)), (19) 

§-^{i\ppim = -i2it){i\ppim. (20) 



In the right-hand side of Eq. (19), the " + " and " — " signs are for bosons and fermions, 
respectively. From the conservation of the number of particles, we have 

7i(t)(l±(/|pp(t)|/))=72(t)(^|Pp(t)K). (21) 

We can set a nonnegative number Wfi(t) so that 71 (t) = Wfi(t){i\pp(t)\i) and 72(0 = 
Wfi{t){l ± {f\pp(t)\f)) in order to satisfy Eq. (21). Then we have 

Af{t) = -lwf,{t){^p,m)\f){fl (22) 

Ai/Ht) = -l^Vf^mi ± (/|Pp(t)|/))|^)(^|. (23) 

Taking A^ (t) and A^f^\t) as Ap{t) and Ap{t), we can derive the master equation for such 
a special case after inserting them into Eqs. (1) and (16). Now Eqs. (19) and (20) become 

^(/Ipp(^)I/) = -|(^Ipp(^)N) = wf,it)mppm)ii ± iflpMf))- 

In the above equation, the change rates for bosons and fermions are different because of 
the last factor, which is determined by the number of particles in |/). For bosons the 
transition becomes faster if there are more particles in |/). On the other hand, for fermions 
the transition is forbidden when |/) is filled with particles, which is consistent with Pauli 
effects. 

In general, the reservoir may induce many transitions when the system contains many 
orbitals. Assume that the transitions are within an orthonormal complete set 5" in which 
each orbital n corresponds to a normalized ket \n). The transition n —^ n' should induce the 
loss and gain in n and n' without changing the number of particles; we shall set A-^'"' (t) = 
-lwn'nit){n\pp{t)\n)\n'){n'\ and A^''''\t) = -iu;„.„(t)(l ± {n'\pp{t)\n'))\n){n\ for such a 
transition. To include all transitions, we shall set Ap{t) = I](n'n) ^p " (^) ^^^ ^p(0 = 
Ein'n) A^/'^Kt)- Eqs. (1) and (16) become 



^^p,it) = -^[H{t),Ppit)]+Y.{Ppit),A^'-\t)}- T.{I±Ppit)^A^'-\t)}, (24) 

(n'n) (n'n) 



= ^[H{t),ppit)] -^Yl «^n'n(0(l ± {n'\Pp{t)\n')){pp{t), \n){n\} 

(n'n) 

+ 2 51 Wn'nit){n\pp{t)\n){I ± pp{t),\n'){n'\}, 

(n'n) 

which is a nonhnear equation. Comparing the above equation to Eq. (A16) derived in 
Appendix A, we can see that the relaxation term can also be obtained formally from the 
second quantization, by considering a system composed of noninteracting particles in contact 
with the reservoir. 



V. DISCUSSIONS 

For stochastic processes, several equations are used under different conditions [1-11]. It 
will be shown in this section that we can derive the Markoff master equation [3-5,7] and 
quasiclassical master equation [9-11] from Eq. (24). An extension of Eq. (24) is also obtained 
after comparing Eq. (24) with the quantum master equation [6]. 

Consider the homogeneous case first. That is, the Hamiltonian and the density matrix 
commute with the momentum operators P, and S is composed of the eigenlevels of P. Let 
|p) be the plane wave with the momentum p. By setting f{p,t) = (p|/9p(t)|p), we obtain 
from Eq. (24) 

p{p, t) = Y: «^pp' (t) [1 ± /(p, t)]/(p', t) (25) 

p' 

-5:«;p,p(t)[l±/(p',t)]/(p,t), 

p' 

which is just the quasiclassical master equation. [2, 10] 

Next we consider the low-density limit. If we diagonalize pp{t) as 

pAt) = J2c,\Mt)){Mt)\, 

X 

then all cx{t) are small, such that / ± pp{t) = J2x[^ ± cx{t)]\ipx{t)) {ipx{t)\ ~ /. In addition, 
A''p-"-\t) in Eq. (24) can be reduced to —^Y.{n'n)'^n'ni't)\n){n\. If coeffecients w„'„(t) are 
independent of t, Eq. (24) is reduced to 



-Ppit) = i[//(t),pp(t)] -^J2 Wn'niPpit), \n){n\} (26) 

(n'n) 

+ Yl Wn'n{n\ppit)\n)\n') {n'\. 

(n'n) 

The above equation is just the so-called Markoff master equation [3-5,7]: 

^Pp(t) = ^[H{t),p,{t)] - 1 ^ «;„'„{Pp(t), \n){n\} (27) 

(n'n) 

+ Ywn'n{n\pp{t)\n)\n') {n'\ - ^ r„/„(n|pp(f)|n')|n)(n'|, 

n'n n^n' 

with all r„/„ = 0. The last term of Eq. (27) is the so-called pure-dephasing term. In the 
Markoff master equation, {r„/„} is a set of nonnegative numbers with r„/„ = Tnn'- It is easy 
to see that for fermions the pure-dephasing term is symmetric with respect to particles and 
holes, and hence it seems that we can add such term into Eq. (1). The positivity of Pp(t), 
however, can be broken under a pure-dephasing term, as shown in Appendix D. 
Eq. (26) is, in fact, a particular case of the quantum master equation [6,12], 

^Pp(t) = ^-[H{t),p,{t)] - i^{pp(t), W,W/} + E W/pp(t)W, (28) 

if we set 

>Vi = v/^Z^|n)(n'|, (29) 

where {wn'n} is a set of positive constants and / = {n',n). Therefore, Eq. (28) can be taken 
as a generalization of Eq. (26). Similarly, we can generalize Eq. (24) to 

^Pp(t) = -^[H{t),p,{t)] - i^{pp(t), WKt)(/±Pp(t))W/(t)} (30) 

+V{/±Pp(t),W/(t)pp(t)WKt)} 
^ I 

with an arbitrary set of operators {VV;(t)}. Eq. (30) can be obtained from Eqs. (1) and (16) 

by setting Ap{t) = -i Ei WKt)(/ ± Pp(t))W/(t) and Ap{t) = -i Ei W/(t)pp(t)WKt), and it 



can be reduced to Eq. (28) in the low-density limit (under which /± Pp{t) ~ J). Therefore, 
Eqs. (1) and (16) are general equations for quantum stochastic processes. It is shown in 
Ref. 13 that under suitable assumptions the positivity of Pp{t) is preserved under Eqs. (24) 
and (30). In addition, for fermions we have {a\pp{t)\a) < 1 for all normalized ket |a), if 
such a relation holds initially. 

As mentioned in section II, we shall include the second term of Eq. (1) if the time- 
evolution operator U{t, t') in Eq. (2) is generated by a nonhermitian Hamiltonian 7i{t). For 
fermions we can operate such a time-evolution operator on the matrix Pp{t), which is defined 
in Eq. (6) for holes, and obtain the following equation 

^p^(t) = ^^[H{t),pp{t)] + {pp{t), Ait)}. (31) 

Replacing A{t) by Ap{t), we can see that the last term of Eq. (1) corresponds to the last 
term of the above equation, after rewriting Eq. (1) as Eq. (13). Therefore, by considering 
the nonhermitian Hamiltonian, Eq. (1) can also be derived by including the last terms 
in both Eqs. (4) and (31), while in section II it is derived from Pauli effects. Setting 
A'{t) = Ap{t) =1= Ap{t) with the signs " — " and " + " for bosons and fermions, actually we 
can rewrite Eqs. (1) and (16) as 

^p,(t) = -^[Hit),ppit)] + {ppit),A'it)} - 2Apit). (32) 

But it is easy to see that for fermions, Eq. (1) is of a symmetric form with respect to particles 
and holes, after rewriting Eq. (1) as Eq. (13). To compare the equations for bosons and 
fermions, we shall use Eq. (16) for bosons, rather than the above one, when Eq. (1) is used 
for fermions. 

In the quasiparticle theory [14-16], which is a good approximation for many-electron 
systems, with suitable assumptions the wavefunction ^(t) of a quasiparticle (or a quasihole) 
is governed by [16] 

^^^mt)) = n{t)Mt)). (33) 

10 



Here the nonhermitian Hamiltonian Ti.{t) is defined in Eq. (3), and its nonhermitian part 
is due to the finite hfetimes of the quasiparticles (quasiholes). The above equation can 
correspond to Eq. (4) or Eq. (31) by rewriting it as 

^p'{t) = -[H{t),p'{t)] + {p'{t),A{t)}, (34) 

at I 

with p'{t) = |^(t))(^(t)|. From the quasiparticle theory, therefore, we can see why the 
nonhermitian part of the Hamiltonian should correspond to the decay of particles or that 
of holes. Since in an orbital the decay of particles (holes) is equivalent to the increase of 
holes (particles), it is natural to consider both the last terms of Eqs. (4) and (31) to obtain 
Eq. (1), as mentioned above. But it should be emphasized that in this paper particles and 
holes are, respectively, the filled and empty parts of any orbitals, while in the quasiparticle 
theory they are the filled and empty parts above and below the Fermi energy. While Eq. 
(1) can be applied to a nonequilibrium Fermi system with no well-defined Fermi energy, in 
the quasiparticle theory the system is not far away from the equilibrium and we shall still 
consider Fermi energy. In addition, Eq. (1) describes the time evolution of a whole system, 
but in the above equation we are only interested in the energy and lifetime of a quasiparticle 
or quasihole, and do not consider the third term of Eq. (1) or Eq. (13). 

VI. CONCLUSION 

It is shown in this paper that the master equation for fermions is of a symmetric form 
with respect to particles and holes, and hence can be derived by extending the mean-field 
Hamiltonian to include an antihermitian part. The equation for bosons is also obtained. 
The derived equations reduce to the Markoff master equation in the low density limit and 
to the quasiclassical master equation for homogeneous systems. 



11 



APPENDIX A: MARKOFF MASTER EQUATION 

Consider a system S coupled to a reservoir R. The total Hamiltonian H = Hg + Hr + V, 
where Hs-, Hji, and V are the operators for the system, the reservoir, and the interaction 
between S and R, repectively. In addition, [Hs,Hji] = 0. If the system is composed of 
noninteracting identical particles, we can write 

n 

with {cn,cl^,} = Sn,n' and {Cn,Cn'} = for femions and [c„,c„;] = 5n,n' and [cn,Cn'] = for 
bosons. 

Assume V = Y^ny^n' /n,n'cJjC„', where /„_„' contain the operators for R and satisfy /„_„/ = 
/„/„. Under suitable assumptions [3,4,7] the density matrix pj^ of the reservoir can be 
expanded by the eigenlevels \r) of H^: 

Pn = J2nEr)\r){r\. {A2) 

r 

Here Er is the eigenvalue corresponding to \r) and F{Er) is the thermal equilibrium distri- 
bution. Let ps{t) be the density matrix of the system and pc{t) be the total density matrix, 
then 

Ps{t) = T.{^\pc{t)\r). {A3) 

r 

In addition, at t = we assume that 

Pc{0) = Ps{0) PR (AA) 

with positive ps{0) satisfying the normalization condition trps{0) = 1. 

Choose 7^0 = Hg + Hf> as the unperturbed Hamiltonian, the time-evolution operator, 
V(^^\t,t'), in the interaction picture satisfies 

^U^'\t,t') = V^'\t)U^'\t,t'), (A5) 



12 



where V^'\t) = En^n' fn,l'{t)e'^'"-'"'^'cicn' with /i5'(t) = e^^«7„,„.e-^^«*. At t > 0, we 
have 

pg)(t)=W(^)(t)pc;(0)W(^)t(^). (^6) 

in the interaction picture. After expanding p^ (t) to the second order with respect to V and 
assume that Trp^ifny = [3,4] we have 

P^s\t) - PsiO) = ^t[psiO), ^r] - t{psiO), S,} + Y. F{Er) f f x (A7) 

where S^ and Sj are the real and imaginary parts of 

E^!^F(E,) fdt, r dt,{r\V^'\ti)\r'){r'\V^'\t,)\r). (A8) 

t^, Jo Jo 

If t is much longer than the correlation time of the reservoir, S is just the self energy [4,8]. 

To calculate the third term on the right hand side of Eq. (A7), let (r|/„^„/|r') = 
\{r\fn,n'\r')\e''"-'^^''"^- Since /„,„, = /^,^„, ^„,„.(r,r') = -^„,,„(r',r). Assume that ^„,„.(r,r') 
is random with respect to r' while |(r|/„^„/|r')| varies slowly with respect to r'. In addition, 
when {n,n') ^ {m',m), (l)n,n',m,m'{i^,r') = ^n,n'(^)^') — ^m,m'(^)^') IS still random. That is, 
the change of On,n'{i^, r') due to the reservoir is random, and there is no correlation between 
different transitions. Under these assumptions, only those terms with (n, n') = (m', m) in 
the summation do not vanish. 

After some calculations, we can set Wn'n = J2r,r' F{Er)\{r\fn,n'\r')\'^SE„+Er,E„,+E^, as the 
transition rate for particles to jump from n to n' so that 

Ignoring T,^ [3,4], Eq. (41) in the Schrodinger picture becomes 
d 1 

g^Ps{t)\t=0 = i[Ps{t),Hs] - 2 II Wn'n{clCn'clcn, ps{t)} (AlO) 

13 



+ Yl Wn'nci,CnPs{t) 



Cj^Cfi' ■ 



Assuming that the above equation holes at any time [3,4] and setting An'n = \/Wn'nCn'Cm 
we can reduce the above equation to be 

which is just the quantum master equation [6] with the second and third terms as the loss 
and gain factors. 

Let \n) be the orbital corresponding to c„, and pp{t) be the one-particle density matrix 
satisfying {n\pp{t)\n') = Trcl^Cn'psif)- From Eq. (AlO), the diagonal terms satisfy 

d 

w;{n\pp{t)\n) = 'Y[-Wn'nTr{clcnCn'cl,ps{t)) + Wnn'Tr{cl,Cn'CnClps{t))], {A12) 

and phases are governed by 

d 

— {n\pp{t)\n") = iiySn" - en){n\pp\n") (A13) 

-9 H '^n'nTr{Cn'c{,ciCn"ps{t)) + o ^ Wnn'Tr{c\,Cn'Cn"c\ps{t)) 



+ 2 ^ Wn"n'Tr{Cn"cic^n,Cn' ps{t)) " o H Wn>n"Tr{c\Cn"Cn>c{, ps{t)), 

where n ^ n" . 



n'^n" n'^n" 



Disjointing from Eq. (A12), so that c„' and c„, are paired, we have 

^Mpp{t)\n) = -Ww,,^il ± {n'\ppit)\n')){n\ppit)\n) (AU) 



+ ^Y.'^nn'{l ± {n\pp{t)\n)){n'\pp{t)\n'), 

n' 

while from Eq. (A13) we have 

— {n\pp{t)\n")\n^n" = t{n\[pp{t),Ho]\n") (A15) 

14 



-^J2(^ri'n" + «'n'n)(l ± {n'\pp{t)\n')){n\pp{t)\n") 



±-^{wnn' + Wn"n') {n\pp{t)\n") {n'\pp{t)\n') . 

n' 

Here the signs " + " and " — " in the symbol " ± " are for bosons and fermions, respectively. 
Although in Eq. (A13) the disjointness is not allowed for n' = n or n" = n, it is suitable to 
ignore such a problem if the summation is over many terms. 

It has been shown that Eq. (A14) corresponds to a semiclassical master euqation. Al- 
though Eq. (A15) looks complicated, with some calculations we can see that pp{t) is governed 
by the following equation, 

O -1 -I 

-Pp{t) = -[H,,pp{t)] - - ^ Wn'n{l ± {n'\pp{t)\n')){pp{t), \n){n\} (A16) 

(n'n) 

+ 2 51 Wn'n{n\pp{t)\n){I±pp{t),\n){n'\], 

(n'n) 

when Eq. (A15) is considered together with Eq. (A14). The last two terms in the above 
equation just correspond to the relaxation term derived in section IV. 

APPENDIX B 

Consider a partially filled band in a crystal. Because holes are vacencies of any orbitals 
in this paper, the reference energy for holes is different from that for particles. In a specific 
orbital, different wavefunctions are used for particles and holes to obtain the correct crystal 
momentum, charge, and the transition energy [10]. We diagonalize Pp{t) as in Eq. (5), the 
number of holes in A is 1 — c\{t), so the density matrix of holes is 

A 

where -0^ is the wavefunction for holes in orbital A. It is easy to see that we can obtain ph{t) 
from Pp{t) after transforming the wavefunctions. 
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APPENDIX C 



To see that we cannot describe the gain of particles in orbital |0) by Eq. (4), we just 
need to consider the special case when no particle occupies |0) at t = 0, i.e. (0|pp(O)|0) = 0. 
For any |0') orthogonal to |0), we have 

mpMW)\ < {<P\Pvm<pf"{<t>'\pMWf' = 0, (ci) 

from Cauchy's inequality, because Pp(0) must be positive. Expanding H{0), A{0), and Pp(0) 
in Eq. (4) by using any orthonormal completet set containing |0), it is easy to see that 
^(0|Pp(O)l0) = for any H(0) and A(0). However, the gain rate should be the largest at 
t = 0. 

The change of the number of particles due to the loss and gain factors should be propor- 
tional to t. If we take those terms of the order t^, of course, the Liouville flow described by 
H[t) can also induce the increase of particles in |0). 

APPENDIX D 

To see that the positivity of Pp{t) can be broken by including the pure-dephasing terms, 
consider a system in which only three orbitals |1), |2), and |3) are allowed, and assume that 
the Hamiltonian Hq is diagonalized by these three orbitals. Let the equation for the time 
evolution of pp{t) be 

^Pp(t) = i[i/o,Pp(t)] -r|2)(2|pp(t)|3)(3| -r|3)(3|pp(t)|2)(2|. (Dl) 

Assume, for example, that initially, 

Pp(0) = ^(|l)(l| + |2)(2| + |3)(3|) {D2) 

10 2 

+-(|1)(2| + |2)(1| + |1)(3| + |3)(1|) + -(|2)(3| + |3)(2|). 

After calculations, it is easy to see that Pp(0) is positive, but one of the eigenvalue of pp 
becomes negative as t — *> cxd. 
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